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Accurate Born-Oppenheimer potential for HeH 
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^ '■ Abstract 

We demonstrate high accuracy calculations for the HeH+ molecule using newly developed analytic for- 
J^ I mulas for two-center two-electron integrals with exponential functions. The Bom-Oppenheimer potential 

for the ground electronic ^S+ state is obtained in the range of 0. 1 - 60 au with precision of about 10^^^ au. 
As an example at the equilibrium distance r = 1.463 283 au the Bom-Oppenheimer potential amounts to 
—2.978 708 310 771(1). Obtained results lay the ground for theoretical predictions in HeH+ with spectro- 
scopic precision. 



PACS numbers: 31.15.ac, 31.50.Bc 



I. INTRODUCTION 

HeH"*" molecule is a simple two electron molecule consisting of the a particle and of the proton 
as nuclei. In the ground electronic ^S+ state, both electrons are mostly centered around a nucleus 
with the proton distance from a being about r = 1.46 au. The first accurate variational calculations 
of the Bom-Oppenheimer (BO) potential of HeH+ reported by Wolniewicz [JJ, Kolos and Peek 
121], Kolos yO, were shortly afterwards refined by Bishop and Cheung 141]. In these calculations 
authors represented the electronic wave function in terms of exponential functions times vari- 
ous polynomials in interparticle distances with general nonlinear parameters (Kolos-Wolniewicz 
functions) or with some parameters fixed (James -Coolidge functions). Accuracy achieved in cal- 
culations by Bishop and Cheung was about 4 ■ 10^^ au. More accurate result ~ 10~® au, at the 
equilibrium distance r = 1.46 au, have been obtained by Cencek et al. pQ using explicitly corre- 
lated Gaussian functions. Even if the whole BO potential is known with their accuracy (~ 0.02 
cm^^), it will be insufficient in comparison with the most accurate (~ 10^^ ^^^l) measurements 
of ro vibrational transition frequencies [|6|]. More recently Adamowicz et al. tu^ performed di- 
rect nonadiabatic calculations for HeH+ including relativistic corrections, however, only for purely 
vibrational states and without estimation of accuracy. 

In this work we use asymptotically correct generalized Heitler-London functions, which are 
the product of an exponential atomic He function times an arbitrary polynomial in interparticle 
distances. Thanks to analytic formulas for the resulting two-center two-electron integrals [lOJ, we 
perform calculations with a large number ~ 20 000 of basis functions. The resulting accuracy of 
about 10^^^ au, or better (for large distances) is 4 orders of magnitude improvement with respect 
to previous values. The accuracy of BO potential is not directly transformed into accuracy of 
rovibrational energy levels. It is because Bom-Oppenheimer energies should be supplemented by 
adiabatic, nonadiabatic, relativistic and QED corrections. As it has been recently demonstrated 
for H2 and isotopomers [11], the perturbative treatment in the electron-nuclei mass ratio up to 
0{m/MY and the fine structure constant a up to 0{aY can provide rovibrational transitions with 
accuracy reaching 0.000 1 cm ^^. Analogous calculations can be performed for HeH+, and results 
obtained here are the first step for spectroscopically accurate theoretical predictions. 



II. ECA BASIS SET 

We represent the ground state nonrelativistic electronic wave function in terms of generalized 
Heitler-London functions, which are the product of the atomic He exponential function with an 
arbitrary polynomial in all interparticle distances, 

= J] q„}(l + Pi2)e-"(^-+'^-VSr;^l<^r^lr,"| (1) 

{n} 

such that 

5 

Y^ni<n, (2) 

for some integer ^2, and call them the explicitly correlated asymptotic (ECA) basis. In the above 
C{„} are linear coefficients, {n} represents a set of 5 numbers {ni, n2, n2, n^, n^}, P12 is the oper- 
ator which replaces ri with r2, indices {A, B) denote nuclei, and (1, 2) denote electrons. Matrix 
elements of the nonrelativistic Hamiltonian can be expressed in terms of integrals of the form 

dVi f d^r2 r e'^'''^ 1 e''^''^^ 1 

,l-"i ».l-"2 r^-n-i l-n4 1-r 
12 ' lA ' IB ' 2A ' 2B 

with nonnegative integers rij. These integrals are calculated as follows. When all nj = the so 
called master integral, which can be derived from a general expression in Ref. [llOll . is 

f{r,f3) 
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where Ei is the exponential integral function. The one-dimensional integral in Eq. dH) is calculated 
numerically using the adapted Gaussian integration. With 120 integration points one achieves 64 
significant digits for all the distances. All /-integrals with higher powers of electron distances can 
be obtained from recursion relations, which were derived in Ref. [llOn . Since 

/(m, n2, ^3, r24, ns; r, (3) = f{n^, ^2, ng, ^4, ^5, 1, /3 r) ^2+n,+n,+n3+n4+n5 (6) 



we present, as an example, formulas at r = 1 for all integrals with ^^ n^ < 2 
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where / = /(I, /3) and F = F{/3). Other /-integrals from the same shell can be obtained by using 
the symmetry /(ni,n2,n3,n4,n5) = /(ni,n4,n5,n2,n3). Integrals with higher powers Ui are 
of analogous form, they are all linear combinations of /, F, Ei, Exp, and identity functions with 
coefficients being polynomials in 1/(3. Using a computer symbolic program, we have generated a 
table of integrals with J2i ^i ^ 37, which corresponds to a maximum value of f] = 16. 

ni. NUMERICAL RESULTS 



The matrix elements of the nonrelativistic Hamiltonian between ECA functions are obtained 



as described by Kolos and Roothan in 1I12I1 . The resulting expression is a linear combination of 



various ECA integrals which are calculated by using analytic formulas as presented in the previous 
section. Their evaluation is fast and accurate, thus suitable for calculations involving large number 
of basis functions. 

Eigenvalues of the Hamiltonian matrix are obtained by inverse iteration method for various 
length of basis sets (Vl = 7 — 16). These eigenvalues, obtained at 132 intemuclear distances 
ranging from 0.1 au to 60 au, form a BO potential. ECA functions are especially suitable for 
large distances, as they include functions for the helium atom alone. As a result, the atomic He 
energy is obtained with an accuracy of about 10^^^ for all the distances, in fact much better than 
that for HeH+. In order to improve numerical accuracy for HeH+, we used a triple basis set 
(n, fi — 2, fi — 4) each one with its own optimized nonlinear parameter a. Numerical calculations 
are performed for r < 12 au using the quadruple precision, and for r > 12 au using the octuple 
precision arithmetics. In order to check numerics, we repeated calculations around the equilibrium 
distance r = 1.463 283 au by using James-Coolidge basis set, for which analytic formulas have 
been developed in Ref. [13]. This basis has a slower numerical convergence but is numerically 
more stable, so we could use even larger number N = 27334 of basis functions. Obtained results 
have similar accuracy and are in perfect agreement with that obtained with ECA functions. 

The most accurate variational result reported at the distance r = 1.46 au is compared in Table 
Uto all the previous results obtained so far in the literature. 

TABLE I: Variational Born-Oppenheimer potential for the HeH+ molecule at r = 1.46 au 
Authors energy [au] 

1965 L. Wolniewicz [IJ -2.978 666 7 

1994 W. Kolos and L. Peek [2] -2.978 689 06 

1979 D.M. Bishop and L.M. Cheung [4] -2.978 702 62 

2006 W. Cencek, J. Komasa and J. Rychlewski [5J -2.978 706 591 

2012 this work -2.978 706 600 341 

In performing extrapolation to a complete basis set (fi — > oo), similarly to the previous H2 
case [13], we observe the exponential e"'^^ convergence, in other words the log of differences in 
energies for subsequent values of f2 fits well to the linear function. This makes extrapolation to 
infinity quite simple. Extrapolated results for the whole BO potential curve in the range 0.1 — 60 
au are presented in Table |ll] 



TABLE II: Numerical values for BO potential at different intemuclear distance r, shifted by the He ground 



state 


energy E = -2.903 724377034120 au. Results 


are obtained by extrapolation to the complete set of 


basis functions, r„^i^ = 1.463 283 au. 
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r/au 
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r/au 
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r/au 


E 


0.10 


15.765 275 621061(3) 


1.60 


-0.072 566 613147(4) 


3.25 


-0.010 505 993 34515(2) 


5.80 


-0.000 684 274 833 49(5) 


0.20 


6.050 585 777 587 9(4) 


1.65 


-0.070 782 049 442(4) 


3.30 


-0.009 805 291916 24(2) 


5.90 


-0.000 634 959 82145(7) 


0.30 


3.035104 004 312 9(1) 


1.70 


-0.068 681986181(4) 


3.35 


-0.009153 886 966 24(2) 


6.00 


-0.000 590184 825 0(1) 


0.40 


1.676 798 340 477 6(1) 


1.75 


-0.066 346 472 331(4) 


3.40 


-0.008 548 577 089 50(2) 


6.20 


-0.000 512 278183 5(1) 


0.50 


0.961224 244 974 8(1) 


1.80 


-0.063 841750 276(3) 


3.45 


-0.007986 304 940 94(2) 


6.40 


-0.000 447 216 535 6(6) 


0.60 


0.551425 9232109(2) 


1.85 


-0.061222 573 648(3) 


3.50 


-0.007464161626 98(2) 


6.60 


-0.000 392 456 5613(6) 


0.70 


0.305137 807 453 6(3) 


1.90 


-0.058 534109 374(3) 


3.55 


-0.006 979 388 833 37(2) 


6.80 


-0.000 346 039 248 5(6) 


0.80 


0.153 061340 049 4(5) 


1.95 


-0.055 813 503 899(3) 


3.60 


-0.006 529 379 042 23(2) 


7.00 


-0.000 306 439 293 9(6) 


0.90 


0.058 077 362 203 3(8) 


2.00 


-0.053 091177182(2) 


3.65 


-0.006111674157 71(2) 


7.20 


-0.000 272 456 6819(6) 


0.95 


0.024 964 448 474 0(9) 


2.05 


-0.050 391895 330(2) 


3.70 


-0.005 723 962 825 47(2) 


7.40 


-0.000 243137 960 03(5) 


1.00 


-0.001080 054 475(1) 


2.10 


-0.047 735 662 718(2) 


3.75 


-0.005 364 076 70123(2) 


7.60 


-0.000 217 718 542 45(4) 


1.05 


-0.021440 483 674(2) 


2.15 


-0.045138 466 546(2) 


3.80 


-0.005 029 985 89214(2) 


7.80 


-0.000195 580 003 55(3) 


1.10 


-0.037 213 297 730(2) 


2.20 


-0.042 612 900 482(1) 


3.85 


-0.004 719 793 768 31(2) 


8.00 


-0.000176 218146 92(3) 


1.15 


-0.049 273 330 916(2) 


2.25 


-0.040168 689 007(1) 


3.90 


-0.004431731314 03(2) 


8.50 


-0.000137 416 959 04(5) 


1.20 


-0.058 323 089 351(2) 


2.30 


-0.037 813130 0178(9) 


3.95 


-0.004164151165 85(2) 


9.00 


-0.000108 799 429 7(1) 


1.25 


-0.064 929 882196(3) 


2.35 


-0.035 551469 970 6(7) 


4.00 


-0.003 915 521460 65(2) 


9.50 


-0.000 087 296 3414(7) 


1.30 


-0.069 554101696(3) 


2.40 


-0.033 387223169 6(6) 


4.10 


-0.003 469 526 010 47(2) 


10.0 


-0.000 070 874 764(5) 


1.35 


-0.072 570 979 391(3) 


2.45 


-0.031322 444 662 9(5) 


4.20 


-0.003 083 563 658 79(2) 


10.5 


-0.000 058152 803(6) 


1.38 


-0.073 739 700 599(4) 


2.50 


-0.029 357 964 4214(4) 


4.30 


-0.002 748 907 76176(2) 


11.0 


-0.000 048170 045(4) 


1.40 


-0.074 287 476 582(4) 


2.55 


-0.027 493 589 050 7(3) 


4.40 


-0.002 458102 55614(2) 


11.5 


-0.000 040 245 989 5(9) 


1.41 


-0.074 497 742 079(2) 


2.60 


-0.025 728 276107 6(2) 


4.50 


-0.002 204 793 867 58(2) 


12.0 


-0.000 033 8901315(2) 


1.42 


-0.074 667 994 434(2) 


2.65 


-0.024 060 285148 3(2) 


4.60 


-0.001983 576 228 38(2) 


13.0 


-0.000 024 539 52511(1) 


1.43 


-0.074 799 949 340(2) 


2.70 


-0.022 487308 864 0(1) 


4.70 


-0.001789 856 934 04(3) 


14.0 


-0.000 018 206 960 42(1) 


1.44 


-0.074 895 257 253(1) 


2.75 


-0.021006 587 043 5(1) 


4.80 


-0.001619 736 634 04(3) 


15.0 


-0.000 013 793 858 65(1) 


1.45 


-0.074 955 506 003(1) 


2.80 


-0.019 615 005 609 25(9) 


4.90 


-0.001469 905 474 75(3) 


16.0 


-0.000 010 641724 04(1) 


1.46 


-0.074 982 223 308(1) 


2.85 


-0.018 309182 582 72(7) 


5.00 


-0.001337 553 485 43(4) 


17.0 


-0.000 008 34146147(1) 


"^min 


-0.074 983 933 737(1) 


2.90 


-0.017 085 542 525 64(6) 


5.10 


-0.001220 293 75156(4) 


18.0 


-0.000 006 630 920 31(1) 


1.47 


-0.074 976 879155(1) 


2.95 


-0.015 940 380 764 26(5) 


5.20 


-0.001116 096 895 45(4) 


19.0 


-0.000 005 337 497 57(1) 


1.48 


-0.074 940 888 099(1) 


3.00 


-0.014 869 918 514 72(4) 


5.30 


-0.001023 235 434 23(4) 


20.0 


-0.000 004 344 79744(1) 


1.49 


-0.074 875 611446(1) 


3.05 


-0.013 870 349 882 66(3) 


5.40 


-0.000 940 236 685 33(4) 


30.0 


-0.000 000 855 679 62(1) 


1.50 


-0.074 782 359 346(1) 


3.10 


-0.012 937881595 48(3) 


5.50 


-0.000 865 843 012 73(5) 


40.0 


-0.000 000 270 476 79(1) 


1.52 


-0.074 516 925 580(3) 


3.15 


-0.012 068 766 236 41(3) 


5.60 


-0.000 798 978 332 42(5) 


50.0 


-0.000 000110 738 51(1) 


1.55 


-0.073 938 987 047(4) 


3.20 


-0.011259 329 676 92(2) 


5.70 


-0.000 738 719 940 4(4) 


60.0 


-0.000 000 053 39146(1) 



IV. SUMMARY 



We have demonstrated applications of analytic formulas for two-center two-electron exponen- 
tial integrals in high precision calculations of Bom-Oppenheimer potential for the ground elec- 
tronic state of HeH+, similarly to previous calculations for H2 1130. The use of EGA basis with as 
much as 20 000 functions provided energies with precision of about 10^^^ — 10^^^ au for intemu- 
clear distances up to 60 au. 

The extension of this approach to calculations of excited states of diatomic molecules such 



as H2 and HeH"*" is not difficult. As long as, the trial function includes at most two different 
nonlinear parameters (in the exponent), the analytic formulas for corresponding integrals are not 
exceedingly large. Such two parameters exponential functions can represent well an arbitrary 
electronic state of a diatomic molecule. In fact, we aim to develop a general code for calculation 
of BO potential of the diatomic molecule, using quad and the arbitrary precision arithmetics. The 
extension to adiabatic, nonadiabatic and relativistic corrections is more problematic. For their 
evaluation one needs integrals with inverse powers of interparticle distances, for which formulas 
can be quite complicated. We think that their derivation for James-Coolidge basis is within the 
reach, therefore for small intemuclear distances, where this basis works quite well, all relevant 
corrections, including yet unknown a^ QED correction can be calculated to high precision, thus 
improving theoretical predictions for H2 and other two-electron diatomic molecules to about 10~^ 



cm ^. 
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